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Effect of weak disorder on delocalization properties of gapped graphene superlattices 
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We study the effect of weak disorder on the delocalization properties of gapped graphene super¬ 
lattice (SL) formed by periodically located rectangular potential barriers. We consider two types of 
the SLs: the SLs with uniform gap and SLs consisting of alternating layers of gapped and gapless 
graphene regions. Using the perturbative approach we obtain an analytical expression for the in¬ 
verse localization length (fLL) derived for the case of randomly fluctuating geometrie and energetie 
parameters, fn the first case, when the barrier (well) width fluctuates around its mean value, the 
corresponding equation for the fLL reveals the presence of the Fabry-Perot resonances, at which 
the localization length diverges. These resonances are exact, i.e., are stored in any degree of dis¬ 
order. It has been found that the localization properties manifest stronger for the particles with 
energies lying in the non-resonant bands where our approach is extremely sensitive to the degree of 
disorder. For the case of weakly fluctuating both barrier and well widths we analytically obtain ILL 
taking correlations into account. The main effect of the correlations, which lead to an increase (or 
decrease) in the localization length, was revealed near the double resonance arising at coincidence 
of two Fabry-Perot resonances associated with barrier and well widths. The random fluctuations 
of the potential strength also lead to the delocalization resonances. However, they exist only in a 
weak-disorder approximation. We found that, for an array composed of alternating strips of gapless 
and gapped graphene modifications these resonances can appear only for normally incident particles 
in contrast to the SL with a uniform gap. For such particles, the delocalization resonances occur 
also in the purely random potential. This means, in particular, that in the one-dimensional case, 
not all the states of the massive Dirac particles are localized in the presence of weak disorder. 


I. INTRODUCTION 


Recent years, much attention of both theoreticians and 
experimentalists has been paid to the graphene-based su¬ 
perlattices (SLs) M. Such interest results from the 
prediction of possible engineering the system band struc¬ 
ture by the periodic potential. This opens different ways 
to fabricate graphene-based electronic devices. An un¬ 
doped graphene is a zero-gap semiconductor. This prop¬ 
erty leads, in particular, to the total transparency of any 
potential barrier for normally incident electrons (an ana¬ 
log of the Klein paradox). At the same time, most elec¬ 
tronic applications are based on the presence of a gap 
between valence and conduction bands. Therefore, it is 
crucial to induce a band gap in Dirac points to control 
the transport of carriers. For this purpose, several ap¬ 
proaches have been studied both theoretically and exper¬ 
imentally. Among them, size quantization in armchair 
nanoribbons, as well as application of external electric 
potentials along the sample edges in zigzag nanoribbons 
were considered [iMl- It has been shown that the gap 
value increases with decreasing the nanoribbon width and 
strongly depends on the detailed structure of the rib¬ 
bon edges. Other proposed mechanisms, which are effec¬ 
tive also in broad graphene sheets, are strain-induced gap 
opening [l3.fl^ , chemical effects of adsorbent atoms and 
molecules [ 1 ^ and substrate-induced band gap forma- 
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tion owing to a breaking the sub-lattice symmetry El- 
The energy spectrum of the Dirac electrons in an epitax¬ 
ially grown on a SiC substrate graphene layer has been 
measured by Zhou et al. 0 . They observed an opened 
up energy gap of about 260 meV in the electronic spec¬ 
trum. It is worth noting also that the k • p Hamiltonian 
of other two-dimensional materials with hexagonal sym¬ 
metry, such as molybdenum disulfide (M 0 S 2 ), is similar 
to the Dirac Hamiltonian for massive particles 0123 . 

Besides, in some publications, various SLs based on 
graphene with spatially inhomogeneous gap (i.e., the par¬ 
ticles mass), and the possibility of their creation are 
discussed |2ll-[^. It was shown, that the spatial mass 
dependence leads to the suppression of Klein tunneling 
and induces confined states [^, . One way of making 

graphene heterostructures with the required gap modu¬ 
lation is a deposition of graphene on an inhomogeneous 
substrate fabricated from different dielectrics. It is also 
possible to use for these purpose an inhomogeneously 
hydrogenated graphene or graphene sheet with nonuni- 
formly deposited CrOs molecules. 

In our previous works [25|, [26| we investigated the 
electronic band structure and transport properties of 
graphene superlattice in which the gap and potential pro¬ 
file are piecewise constant functions. It was shown that 
in such SL, up to some critical value of potential Vc , 
allowed subbands are separated by gaps. At U > Uc the 
contact or cone-like Dirac points appear in the spectrum. 
It was also found that each a new Dirac point manifests 
itself as a conductivity resonance and a narrow dip in the 
Fano factor F similarly to a gapless SL. However between 
the resonances, behavior of the Fano factor in the con- 
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sidered structure is more complicated and differs from 
pseudo-diffusive behavior (F = 1/3) typical for a gapless 

SL 

Meanwhile, real graphene super lattices cannot be per¬ 
fectly periodic due to random imperfections resulting, for 
example, from variations of the system parameters such 
as potential height, gap value, potential width or barrier 
spacing. It is well known, that in the presence of white- 
noise disorder all the electronic states are localized in the 
thermodynamic limit for a traditional semiconductor su¬ 
perlattice. On the contrary, a sample of gapless graphene 
in the presence of a random one-dimensional potential be¬ 
comes completely transparent for the normally incident 
particles regardless of the sample length and strength 
of disorder. This means, that the states of the mass¬ 
less Dirac particles are entirely delocalized for arbit rary 
disorder strength due to the chiral symmetry [^, 1^ . 
The transport properties of disordered graphene super¬ 
lattices have been studied by several groups jsol-ls^. It 
was found that the transport and spectral properties of 
gapless graphene superlattices created by applying ei¬ 
ther periodic or disordered smooth scalar potentials are 
strongly anisotropic. The dc conductance of graphene su¬ 
perlattice consisting of p-n junctions for various strengths 
of structural disorder imposed on the material has been 
investigated numerically in Ref. [^. It was shown that 
there exists a range of angles around the normal in¬ 
cidence angle, for which the transmission becomes fi¬ 
nite in the presence of structural white-noise disorder. 
For weakly disordered both scalar-potential and vector- 
potential graphene SLs the localization behavior of mass¬ 
less Dirac particles was studied in Ref. [32[ numerically 
as well as analytically by a weak-disorder expansion. In 
particular, strong dependence of the Lyapunov exponent 
(the inverse localization length) on the incident angle of 
the charge carriers injected to a graphene superlattice has 
been predicted. The effects of gap fluctuations on trans¬ 
mission and conductance of the monolayer and bilayer 
graphene SLs were treated numerically in Ref. [33| . 

The aim of this work is to study the effect of weak dis¬ 
order on the localization length and transport properties 
of disordered gapped graphene SLs including the sam¬ 
ples with spatially inhomogeneous gap. We extend the 
theoretical study developed earlier for the case of gapless 
graphene SLs [33 and obtain an analytical expression for 
the localization length derived for the cases of randomly 
fluctuating parameters of the SLs. This expression is in 
a good agreement with direct simulations. We also take 
into account possible correlations for the case of weakly 
fluctuating widths of layers forming the unit cell of the 
superlattice. 

The paper is organized as follows. Section II is de¬ 
voted to the description of the model and the method. 
In Sec. Ill we present the dispersion relation and trans¬ 
mission for graphene-based multibarrier periodic struc¬ 
ture with spatially dependent gap in the presence of the 
step-like potential. The analytical expression for the 
inverse localization length of periodic-on-average disor¬ 



FIG. I: Models of a periodic structure: A = Ao - for the ho¬ 
mogeneous superlattice (HSL); Ao = 0 - for the superlattice 
formed by the alternating layers of gapless graphene (MSL). 


dered graphene SLs as well as the results of numerical 
simulations are presented in Sec. IV. We make a sum¬ 
mary and concluding remarks in Sec. V. 


II. MODEL AND METHOD 

We consider the propagation of an electron through the 
lateral disordered structure formed by a sequence of N 
barrier regions with width dn{Ti = 1,2,... N) separated 
by inter-barrir distance (wells), as shown in Fig. [H 
The disorder is introduced as random, small variations 
of the barrier strengths or other barrier characteristics 
(e.g., gap value in gapped graphene) as well as the barrier 
and well thicknesses around their mean values. In the 
absence of the disorder the considered system is periodic 
with the period / = a + d, with d = (d^), a = (an)- The 
main subject of our study is the localization length Lioc, 
defined as 


7 = 


Lio 


lim 

N^oo 


InTjy X 

2N /’ 


( 1 ) 


where T/v is the random tramsmission coefficient of a 
sample of the length Nl and angular brackets are used to 
denote averaging over different disorder realizations. To 
calculate the transmission coefficient, we use the common 
transfer-matrix approach. 

In the general case, between barriars, where V{x) = 0, 
the wave function 2 pri{x,y) can be written as = 

e*^y^'0n(^), where is a superposition of the wave 

functions of right- and left-moving particles 

tpnix) = Xn + dn < X < a;„+i(2) 


and V are spinor amplitudes, defined by the specific 
Hamiltonian. The amplitudes 7 /+ = and 7 /“ = 

in neighboring wells are mapped from n to 
n + I by the transfer matrix for a single unit 


V^n+l 


= Mn 



(3) 
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with 


Mn 


Here reflection (r^) and transmission (t^) amplitudes are 
determined by the parameters of the nth barrier, as well 
as dynamic characteristics of the particles, i.e., energy 
E and momentum hky. By construction, the transfer 
matrix across N barriers is the product 

N 

Piv = n (5) 

n=l 



so the transmission probability 


Tn 



(6) 


For infinite periodic structure with period I = a d 
matrix elements Mn does not depend on n and the elec¬ 
tronic band structure of the corresponding superlattice 
is governed by the following relation 


2 cos Kl = TrM 


(7) 


with K the Bloch wave vector. As known (see, e.g.. 
Ref. [l^), the transmission coefficient of the array 
of N identical cells (Fig.[T]) can be written in closed form 
as follows 


rpO _ 

J-N — 



rp 2 

t 


sin Nt] 
sin 7^ 



(8) 


where the Bloch phase rj = Kl is defined by Eq. O- 
Note, that a similar expression holds for the transmis¬ 
sion coefficient through an array of length L = Nl 
bounded by the regions, characterized by barrier param¬ 
eters. In this case T'% depends on the reflection (r') 
and transmission {t') amplitudes across the inter-harrier 
region (well). 

In the presence of weak disorder the localization length 
can be found by substituting Eqs. ® ® into definition 

m and expanding the logarithm up to quadratic terms 
in disorder. To perform these calculations, we need to 
know the transmission and reflection amplitudes for a 
single barrier. 


III. GAPPED GRAPHENE SUPERLATTICES 


System under consideration is graphene-based multi¬ 
barrier structure with spatially dependent gap A{x) in 
the presence of the step-like potential V{x). Eor a peri¬ 
odic system V{x) = U(x + /), A{x) = A{x + /), where 


V (x), A{x) 


U, A for 0 < X < d, 
0, Ao for d < X < 1. 


(9) 


The wave functions of charged particles in such model 
obey the Dirac equation with the Hamiltonian 

H = VFPcr A{x)az E V{x)^ (10) 

with p - the momentum operator, a = (<Ja,,cr^), az the 
Pauli matrices and vf ^ 10^ ms“^ the Eermi velocity. 
The equation {H — E)ip = 0 admits the plain wave solu¬ 
tion of the form %p{x^y) = with 


( Q+iky ) 

( -qXiky ^ ■ ( 11 ) 

Eor convenience, hereafter all lengths will be expressed 
in the units of the mean period 1. The natural energy 
scale is Esl = {Esl ~ 6.25 meV for I = 100 nm), 
so that other dimensionless parameters are e = EjEsLi 
v{x) = V{x)/Esl, A(x) = A(x)/Esl, where V(x) and 
A(x) in the barrier and well regions are determined by 
Eq. (j9j). Q is the dimensionless wave vector along the 
x-axis 


Q = 



for barrier, 
for well. 


( 12 ) 


where 


= ^(e - v)^ - A2 - kl, kx = - Ag - fcg. (13) 

By applying the continuity of the wave function at the 
boundaries, we obtain the transmission (t) and reflection 
(r) amplitudes for electrons incident at an angle do with 
respect to the x-axis (Eig. [T]) 


1 

t 


= cos yd + 7- 


£V + AAo - £^ + kl 

—- sin/3, 


kxQx 


(14) 


t kkxQx 

Ei (Aq - AAo - sv) ky 

kkxQx 


(15) 


where fd = Qxd^ k = “ ^o particle wave vec¬ 

tor outside the barrier, do = tan“^ Then tunneling 
through a single harrier is given by 

T (e, ky) = E = (1 + if (e, ky) - l) sm' (16) 


with 


/ (e> ky) 


ev + AAo — e'^ + ky 
kxQx 


(17) 


Transmittance through a single well differs from this ex¬ 
pression only by replacing ^ ^ a = kxa. AtA = Ao = 0 
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the expression (US coincides with the similar to gapless 
graphene [35|, and for Ap = 0 , A 7 ^ 0 is the same as 
that established in Ref. [^. Using Eq. (jS]), we obtain 
the transmission across N identical barriers 


{s, ky) 


^1 + {f {e, ky) - 1 ) sin^ /3 


/ sin Nri\ 
\ sin 7 ^ ) 


(18) 


where Bloch phase 77 according to Eqs. o, o, da can 
be obtained from dispersion relation 


cos 77 = cos a cos ^ f {e^ky) sin a sin /3. (19) 

As follows from the expression m the transmission 
T^{e,ky) = 1 for any N under conditions sin/3 = 0, 
/3 7 ^ 0 or when sin A 77 / sin 77 = 0. The first equation 
Qxd = TTiJi determines the Eabry-Perot resonances & 
related with the barrier regions and the second pro¬ 
duces N — 1 Eabry-Perot oscillations in each allowed en¬ 
ergy band. Similarly, for particles incident on array of N 
unit cells from the barrier region T'% ( 5 , ky) = 1 when 
sin a = 0 , (a 7 ^ 0 or sin Nr]/ sin 77 = 0. Note also, that 
for gapless superlattices with A = Ap = 0 , gapped with 
a uniform gap (homogeneous or HSL) with A = Ap 7 ^ 0 
and for the super lattices, formed by alternating strips 
of gapless and gapped graphenes (mixed or MSL) with 
Ap = 0 , A 7 ^ 0 the function f (e^ky) (pTj) determining 
transport properties and the energy spectrum of the SLs 
has the same form 


to the assumption. Numerically, for generating random 
sequences we use the flat distribution on a finite inter¬ 
val [—^, 6]. An analytical expression for the inverse local¬ 
ization length (ILL) 7 can be obtained by the method of 
perturbation theory. To do this, follow Zhao et.al. [32|, 
represent the expression for the transfer matrix for a sin¬ 
gle unit dH) as 


/ ^iPn \ 


(23) 


where sin(p^ = jr^j, cos cpn = \tn \ and the parameters rrin 
and Pn are determined by the value of kxa and the phases 
of the reflected and transmitted waves (Eq. (j4])). Then 
the weak-disorder ILL or Lyapunov exponent depends 
only on the parameters of the underlying regular array 


7^ = 


5^(7? 


■ tan^ if 


,2 ( siuTn 

P + . 

simp 


tan^ (p 
sin^ 77 


(24) 


where the prime (... )^ denotes differentiation with re¬ 
spect to the perturbation variable s. Using Eqs. (HI, 
m it is not difficult to show that for the considered 
gapped SLs 

tan^ p = (/^ ( 5 , ky) — 1 ) sin^ /3, (25) 


sinm sin a cot /3 — f {s, ky) cos a 

“ 7/2 {e,ky) - 1 


(26) 


p = a---0o-^, 


(27) 


/(£,fc,) = ^^, (20) 

where wave vectors kx and qx are defined by formula (USD 
for each type of the superlattice. 


IV. LOCALIZATION LENGTH FOR 
DISORDERED GRAPHENE STRUCTURES 

In what follows, we consider the disordered multibar¬ 
rier graphene structures in which the disorder is caused 
by random fluctuations of barrier strength, or gap magni¬ 
tude (inside the barriers) as well as by random variations 
of both barrier and well widths. Specifically, we assume 
a weakness of both types of disorder 


^ = tan ^ 


/ Ag - AAq - ev 

V^(A — Ao) + 7 ;Ao 


tan ^0 


(28) 


Note that Eq. ([24j) is correct inside the energy bands 
(TrM < 2) apart from the band edges 77 = 0 , 7 r. When 
TrM > 2, the energy lies in the forbidden miniband. 
In this case ILL 7 ^ is defined by A+, the largest of two 
eigenvalues of the transfer matrix M 7 ^ = In | A+ 1 . 

Using the weak-disorder approach it is possible to gen¬ 
eralize the expression for ILL (|M|) to the case when the 
correlations of fluctuating quantities (e.g., the geometric 
parameters of the structure) exist. Details of the calcu¬ 
lation are given in the Appendix. 


A. Geometric disorder 


Sn = S (1 -1- p^), 5 = 7 ;, A or d, a. ( 21 ) 

Here the index n enumerates the nth unit (d, a) cell, 
is random uneorrelated variables with zero average and 
small variances (j‘1 ^ 1 , i.e. 

(dn) ~ i^PnPn'^ ~ s^nn'^ss' ‘ (22) 

The averaging (...) is performed over the whole array of 
layers or due to the ensemble averaging, that is equivalent 


Eor weakly fluctuating widths of layers (positional or 
geometric disorder), according to Eqs. (|24|) . (jA9p the Lya¬ 
punov exponent, which includes the correlation term, can 
be written as 


7a,d = {f (s^Oo) - 1) 


[ct^CTa P ^ ^ 

2 sin^ T] 


—2al3aad sin a sin f3 cos p] 
2 sin^ p 


( 29 ) 
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FIG. 2: Inverse localization length 7a versus particle energy £ at /cy = 0 for graphene MSL (a, b) and HSL (c, d) corresponding 
to the fluctuation of distance between the barrier with two disorder strength: S = 0.05 (a), (c) and S = 0.2 (b), (d). For all 
cases a = d = 0.5, 'U = tt, A = tt/S. The continuous (blue) curve corresponds to the analytical results [Eqs. ([32]), (l33]) ]. dotted 
line presents numerical data for the structure length A = 10^ with an ensemble averaging performed over 80 realizations of 
disorder. The continuous (red) line shows the dependence of ILL on s in the forbidden minibands, shaded in Figure. Vertical 
lines mark the delocalization resonance positions, placed in the resonance zone shown in the insets. 


where f (e^ky) is given by Eq. (|2Q|) and two types of the 
considered gapped graphene structures differ only in the 
value of the wave vector in well region kx = k^^^ = 

^— /S? — ky for the SL with uniform gap and kx = 
~ with piecewise constant 

gap. 

Assume that only the distances between the barriers 
display random fluctuations around their mean value a, 
that is crj = (Tad = 0, cTfl ~ and 5 determines the de¬ 
gree of disorder. In this case Eq. (1211) indicates that the 
localization length turns into infinity when performing 
the Fabry-Perot resonance conditions Qxd = irm. More¬ 
over, this result is exact, that is valid for any degree 
of disorder. Indeed, in this case transmission amplitude 
across a single barrier t = 1 dH resulting in total trans¬ 
parency of the array from N identical barriers located 
randomly. Figure [2] compares the analytical results for 
7 ^*^^ and 7^*^^ (left and right column, respectively) 
with the numerical data, for a varying particle energy 
5, fixed 6>o = 0 and for two different values of disor¬ 
der: ^ = 0.05 (Figs. |2fa), |2fc)) and 5 = 0.2 (Figs. |2]^b). 


[2](d)). The potential strength is = tt and the gap value 
is A = 7r/3. Numerical data are obtained for an ar¬ 
ray size A" = 10^ with additional average over 80 real¬ 
ization. From Figs. Ha), |2]^c) it can be seen that the 
theoretical expressions for the inverse localization length 
for both models of the SL provide a very good descrip¬ 
tion for the case of weak disorder {6 = 0.05) for this 
range of energies. For higher disorder strength {5 = 0.2) 
the agreement between the analytical and numerical cal¬ 
culations remains good for the SL with a uniform gap 
at least for energies 5 > 8 (Fig. Efd)). But the lattice 
formed by alternating stripes of the gapped and gapless 
graphene (MSL), is more sensitive to fluctuations of the 
inter-barrier distance (Fig. Efb)). It is also clearly seen 
that in the neighborhood of the Fabry-Perot resonances 
the localization is strongly suppressed and the greater 
the amount of the disorder, the narrower this neighbor¬ 
hood. These results are in complete correspondence with 
those obtained (and has been confirmed experimentally) 
in Ref. [39| , which deals with the propagation of electro¬ 
magnetic waves through one-dimensional disordered bi¬ 
layer structures whose unit cell consists of two different 
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FIG. 3: Inverse localization length 7a,d versus particle energy 
£ for uncorrelated (a), completely correlated (b), and anticor¬ 
related (c) disorder at = cos“^ 2/3 for graphene MSL with 
= Stt, A = 27r corresponding to the fluctuations both the 
distance between the barriers and their width near an aver¬ 
age value a = d = 0.5 with disorder strength 6 = 0.005. The 
solid line represents the analytical results, dotted line corre¬ 
sponds to the numerical simulations for an array of A = 5T0^ 
(iV = 3 • 10^ for insert) with the average over 100 realizations 
of disorder. The vertical solid lines indicate the positions of 
the Fabry-Perot resonances arising when qxd = 7r;37r. The 
vertical dashed lines depict the positions of the double res¬ 
onances that occur when Qxd = 27r and fca = tt. Insert 
illustrates the behavior of 7^/^ in a neighborhood of double 
resonance. 


dielectrics. Applicability approximate expressions (ESI) at 
a given disorder strength depends on of the lattice param¬ 
eters. Thus, with increasing the gap magnitude visible 
discrepancies with the results of numerical calculations 
are observed even at the disorder strength (5 = 0.05. 

Similarly, for fluctuating barrier widths (cr^ = a ad = 0, 
crj = /3) the expression ([29|) for the inverse localization 

length ^d reveals the delocalization resonances for both 
types of the SLs under the conditions sin a = 0. When 
these conditions are met, the transmittance T'at through 
N identical wells, separated by barriers, the widths of 
which vary randomly, is equal to one for any N. On the 
other hand, the transmission coefficient T/v through N 
regularly spaced barriers with fluctuating width does not 
decrease exponentially with increasing N, which leads to 
7 c^ = 0, i.e. the suppression of localization. 

For the array with randomly varying both barrier 
width and inter-barrier spacing the ILL ^a,d is obtained 
from the general expression (|2^ . As above, we take 
a = d. Then for uniform random perturbation with the 
same amplitudes on both layers, we have = crj = (5^/3. 
If the disorder is uncorrelated, then a ad = 0- For com¬ 
pletely correlated disorder, when the barrier and well 
widths in the period change in the same way, we have 
(^ad = S‘^/3. In the case of completely anticorrelated 
disorder period of the SL remains constant and we take 
O'ad = —S‘^/3 [i^. For the Fabry-Perot resonances occur¬ 
ring when a = irn OT P = irk with n,k = 1; 2; 3;..., the 


factor sin a or sin /3 in Eq. (E3) vanishes, so that the cor¬ 
relations do not affect the localization properties of the 
structure. It may happen that some of the resonances 
due to different graphene layers coincide for certain val¬ 
ues of 5 and Oq that it is possible under the condition 
sin a = 0, sin (3 = 0. But under these conditions the de¬ 
nominator ^a,d vanishes also m, so localization length 
remains finite and its value significantly depends on the 
existing correlations. In Fig. [3] the ILL for: (a) uncor¬ 
related, (b) completely correlated and (c) anticorrelated 
disorder and oblique incidence is shown as a function of 
energy for disorder strength S = 0.005 for graphene MSL. 
In this Figure the central band includes the value of the 
energy cq = 37r which determine the position of the dou¬ 
ble resonance defined by the conditions a = tt, /3 = 27r. 
Fabry-Perot resonances corresponding to the conditions 
P = TT and (3 = 371, for the chosen lattice parameters 
occur in the first and third allowed energy bands. Anal¬ 
ysis ([29|) shows that the values of 7^/^ corresponding to 
such resonances close to the minimum values of the ILL 
in the relevant energy bands, but do not coincide with 
them. It is interesting also to note that for the SLs with 
spatially inhomogeneous gap the resonance values of the 
ILL does not depend on the energy 




(A^ + v‘^ tan^ So) 
2 


(30) 


Fig. [3] clearly demonstrates a significant discrepancy of 
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FIG. 4: Inverse localization length 7 ^; as a function of energy e for graphene MSL for an angle of incidence = 0 (a, b) and 
HSL for ^0 = 7r/6 (c, d) with fluctuating barrier height for two disorder strength: S = 0.01 (a, c); 5 = 0.5 (b, d). For both cases 
a = d = 1/2, 'L) = TT, A = 7r/3, the other notations are the same as in Figured Numerical calculations were performed for an 
array composed of A = 5 • 10^ until cells {N = 2 • 10^ for inserts). 


approximate analytical calculations with numerical data 
in the first and third energy bands, even with such a rela¬ 
tively small amount of disorder, as (5 = 0.005. In the cen¬ 
tral energy band, containing double resonance, our ana¬ 
lytical predictions, based on the formula ([29]), are more 
consistent with the numerical data. Clearly seen, that the 
most influence of correlations appear close to the double 
resonance: completely correlated disorder suppresses the 
localization length (Fig.ISfb)) and anticorrelated disorder 
term in Eq. (|29]) leads to its increasing (Fig. [3l^c)). For 
the superlatice with uniform gap the results are similar. 


expressions 


MSL ^ ^ 

ky/S? sin^ P [F^ sin a + cos a)" 


klql^m T] 


m 


where 


S = v‘^ky -f- A^kl, 


(32) 


pM 


d{v 



V — e 


Qx 


sin P cos /3, 


(33) 


B. Compositional disorder 


GM 



vkl {ev - kl) \ 
^ ) 



(34) 


.,HSL 


v^a^Z 

sjn2 


(F 


H 


sma - 




COS (a 


In the case when the barrier height fluctuates around 
their mean value the Lyapunov exponent 7 given by 
Eq. (|M|) for both types of superlattices is defined by the 


and 


2 


(35) 
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with 


Z=(A2+fc2), (36) 


= v{v — e)d — ——— sin [3 cos /3, 

Qx 

G" = sin^ /3. 


(37) 


Expression (ESI) takes the simple form in the limiting case 
of (^-function barriers (that is when d ^ 0, v ^ oo but 
such that their product remains constant vd = ip) 


with 


Jcp 


d‘^Z sin^ a 
2kl sin^ T] ’ 


5 . 

cos T] = COS a cos (/:? + — sm a sm ip. 

kx 


(38) 


(39) 


Here ip denotes the mean value of the “potential” at the 
n-th site: ipn = ip F Sipn, where 6ipn are homogeneous 
random perturbations and = {Sip‘^). 

In this case, the dependence of ILL on the barrier char¬ 
acteristics {v and d) are determined only by means pa¬ 
rameter ip therefore, the formula (|38|) can also be ob¬ 
tained as a limit of the expression (ESI) provided that 
only barrier width fluctuates. It is obvious that in the 
limit of very narrow {d 0) barriers superlattice with 
non-uniform gap becomes gapless and the equation (|38|) 
at A = 0 coincides with that of Ref. . The presence 
of a gap leads to localization of the particles incident on 
the structure at arbitrary angles, with the exception of 
delocalization resonances (a = 7rn, n = 1, 2,...) that as 
in the case of gapless SL are exact for arbitrary disorder 
strength. When ip = 0, that is for purely random S- 

potential Eq. EH) reduces to 7^*^^ = , which means 

complete localization of massive Dirac particles in such 
structure. 

Now we return to the general case of a rectangu¬ 
lar potential super lattices. As can be seen from the 
equation m, for array with non-uniform gap all the 
states with ky ^ 0 (i.e., in the case of oblique inci¬ 
dence of the particles) are localized. When incidence 
angle = 0 perturbative delocalization resonances are 
determined vanishing expression in round brackets in 
Eq. (|3T|) . but unlike the geometric disorder, they ex¬ 
ist only for weak disorder strength (Eigs. SJ^a), SJ^b)). 
The weak-disorder expansion for the Lyapunov expo¬ 
nent for graphene SLs with uniform gap Ei) also 

manifests an emergence of the delocalization resonances 
(at any angles of incidence) that occur under condition 
sin a + cos a = 0 and disappear with increasing 

disorder strength (Eigs. lU^c), lU^d)). Eig. |4] also clear 
demonstrates that although with the growth of disor¬ 
der resonances disappear, but near them approximate 
formulas (j35j) describes the ILL well even at a high 
degree of disorder S = 0.5. Note that for a random poten¬ 
tial without a regular superlattice component (vn = Sv^) 



FIG. 5: Inverse localization length 7 ^'^^ at ^0 = 0, A = tt/S, 
a = 1/3 and d — 2l?> versus energy e for a purely random 
potential {v — 0 ) with the amplitude of fluctuations 5v — 
tt/IO. The analytical results [Eq. (00])] and numerical data 
are shown by, respectively, the solid (blue) and dotted curves. 
Vertical lines mark the resonant energies. 


and for ky = ^ the Lyapunov exponent Eq. E3 reduced 
to 


.,HSL 


72 sin2 j 3 




= (,5.2), 


(40) 


Unlike the case of (5-function barriers, this expression 
vanishes for the resonance energies 


e„ = ±yA2 + (^)/ n = l,2,... (41) 

leading to divergence of the localization length (Eig. [5|). 
This contradicts the statement of Ref. j29[ that, in one¬ 
dimensional case massive Dirac particles should be local¬ 
ized for any weak disorders. 

In the case where a fluctuating parameter is the gap 
value in the barrier region (and in the intervals between 
the barriers still A = 0), we find from Eqs. (|M|) - (|28|) 


^MSL 


sin^ P 


f /S? sin^ {3 
I klql sin2 7 ? 


(Fa sin a + Ga cos a)^ + v‘^ky 



Fa 


PS — {ev — /c^)^ sin/3 cos p 

kxQx sin^ P 


(43) 


Ga = £v — k^. (44) 

Similarly to the previous case of fluctuating barrier 
height, the weak-disorder resonances are possible only 
for normally incident particles. Note that although these 
resonances obtained in the weak-disorder approximation 
they survive at significant deviations of the gap from its 
average value. Thus, numerical simulations carried out 
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for the MSL with v = it^ A = 47 r/ 9 at/c^ = 0 showed 
that at 20 % of the gap fluctuation resonance remains well 
defined, although in the non-resonant zones there is a sig¬ 
nificant discrepancy between the approximate analytical 
( 02 ]) and numerical results. 


V. CONCLUSION 

In summary, we have studied the localization behavior 
of Dirac particles in disordered graphene super lattices. 
Using the weak-disorder approximation, we obtained the 
analytical expression for the inverse localization length 
(Lyapunov exponent). The main attention was drawn to 
the two models of the SLs. One of them corresponds to 
the massive Dirac particles (the SL with homogeneous 
gap) in the presence of one-dimensional piecewise con¬ 
stant potentials. Another discussed model is a layered 
structure, made of gapped and gapless graphene strips. 
It is obvious, that the presence of a gapped graphene frac¬ 
tion in disordered SLs leads to suppression of the Klein 
tunneling and localization of Dirac particles with zero 
incidence angles. When the disorder emerges due to ran¬ 
dom thickness variations in the well (or barrier) layer the 
Fabry-Perot resonances leading to divergence of the lo¬ 
calization length arise in the discussed SLs like in other 
models of the Kronig-Penney type. This result holds for 
the SLs, i.e. for infinite systems. When we consider a 
transmission probability through the lateral structure of 
finite length L = Nl, composed of alternating barrier and 
well strips (with the same or different gap magnitudes), 
the outer regions (x < 0 and x > L) may correspond to 
both the well and the barrier parameters. Obviously, the 
transmittance will depend on these boundary conditions: 
T/v - for the well outer regions and T'at for the barrier 
ones. For example, when the width of wells fluctuate 
and resonance conditions Qxd = irn are fulfilled, each of 
the barriers (and hence N barriers) becomes completely 
transparent, i.e. Tat = 1 for any N and for any strength 
of disorder. At the same time, the equation 7 = 0 means 
that corresponding transmission probability through N 
wells T'tv cannot be an exponentially decaying function 
of N (numerical data indicate that T'n does not decrease 
for a system made up of a sufficiently large number of lay¬ 
ers). We also received the analytical expression for the 
localization length for the case of weakly fluctuating both 
barrier and well widths, taking into account possible cor¬ 
relations in disorder. We have studied and compared the 
cases where disorder is uncorrelated to cases where it is 
entirely correlated and anticorrelated. The main effect of 
the correlations, leading to an increase (or decrease) in 
the localization length, was found in the vicinities of dou¬ 
ble resonance arising under the conditions sin kxCi = 0 , 
sinqxd = 0 . 

Also, delocalization resonances for both types of the 
SL are obtained for the barriers with randomly varying 
height, but in contrast to the Fabry-Perot resonances, 
they are approximate. Resonance values of energy and 


angle of incidence are determined by the parameters of 
the system and, in general, can be found only numer¬ 
ically. Corresponding expressions for 7 ^ im. (|3SD) 
demonstrate distinct features of two superlattice mod¬ 
els: for massive Dirac particle resonance condition can be 
performed at arbitrary angles of incidence, while in the 
structures with non-uniform gap such weak-disorder delo¬ 
calization is possible only for the Dirac particles with zero 
incidence angle. Interestingly, the delocalization states 
exist in one dimension (i.e. when ^0 = 0) for the massive 
Dirac particles with energies Sn (|4T]) placed in the purely 
random potential i.e., with the barrier height Vn being a 
constant, randomly distributed in a certain range, which 
determines the degree of disorder. At the same time, 
for the disordered (5-function potential without a regular 
superlattice component all the states are localized. 

The obtained results for the localization length can 
be used for finding an analytical expression for the 
transmission coefficient through the finite-size disordered 
graphene system. In turn, it is possible to apply this ex¬ 
pression to an analysis of conductance and resistance. 
However the analytical results were obtained by statis¬ 
tical averaging. Therefore, their comparison with cor¬ 
responding experimental data requires sufficiently long 
structures. Despite this limitation, such a comparison 
for a single-mode microwave waveguide composed of rel¬ 
atively small number of cells {N = 26) has shown quite 
good agreement in the case of weak disorder for the 
whole range of frequencies . In graphene heterostruc¬ 
tures, a finiteness (or rather smallness) of mean free path 
and phase coherence length is another limiting factor 
for the emergence of sizable disordered superlattice ef¬ 
fects (in particular the delocalization Fabry-Perot reso¬ 
nances). Accordingly, samples with high carriers mobil¬ 
ity and small superlattice period are needed. 
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Appendix A 

To study the effect of disorder correlations we allow the 
width of barriers as well as the distance between them 
fluctuate relative to their average values, adding to the 
relations ([ 22 ]) correlator {PnPn')’ 

77) =0, {Ptpt') = CrjSnn', 

{PnPn') ~ ^a^nn’i {pnPn'} ~ ^adorin'- 


(Al) 
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Accordingly the weak-disorder expansion of the transfer- 
matrix Mn is 

Mn = M + aM'^pl + pM'ppi + yML {plf 

+ {pi)^ + apM::ppyi. (A2) 

Here the matrix M 0 and derivatives of M with re¬ 
spect to the fluctuating super lattice dimensions (M' and 
M") determined by the parameters of unperturbed SL. 
Following Ref. [32|, we find the matrix element (P/v)ii 5 
using the representation (M-representation) in which the 

matrix M is diagonal M ^ M =diag(A+, A_) 


Here the terms, proportional to and lead to the ex¬ 
pressions for 7a and defined above by Eq. (|2^ without 
taking into account the correlation term a/3). To find 

it, we need to know and Computing respec¬ 

tive derivatives of the matrix M and making the trans¬ 
formation to the M-representation, after some algebraic 
calculations, we have 



A+ sinr? ’ 

e-^^'Re - b'*b'pj + Rea^ 

i sin 7 


(A5) 

(A6) 




N 


= A 


N . 


Aj 


Y,\aM'„pl + P^pt 


k=l 


+ ~^^aa (Pk)^ + {Pk) 


where a and b denote matrix elements of matrix M 
(Eq. (jl])): a = Mu, b = M 12 . Using Eqs. ([T4|) . ([15]), 
we get from (|A5p and (|A6p the expressions that define 
the correlation term 


-^al3M'J^/dp%pi 


(A3) 


Note that in this expression, we have omitted the terms 
proportional pkpk'^k^ that vanish in the subsequent av¬ 
eraging. In order to evaluate the disorder-induced Lya¬ 
punov exponent one need to combine Eqs. (HD, (0 and 
(jASp and expand the logarithm within the quadratic ap¬ 
proximation in the perturbation parameters. Performing 
averaging with help of Eq. m, we obtain 



(sin a cos /3 
sin^ 77 


—/ cos a sin P) (— sin P cos a-\- f cos P sin a) 


sin^ p 


(A7) 


7= 


(*;) 


11 


“h Re 




-2Re^ 


(«i) 


TRe 


(«;) 


11 


(*«) 


("«) 


A+ 

2 


11 


2Re" 


11 


^2aPaad 


dm 


(*«) 


Re- 


11 




(A4) 



(A8) 


Next, substituting Eqs. dMD, (EHD into (1X41) we obtain 
the contribution of the correlation term A^cor Ihe in¬ 
verse localization length 


Aye, 


(l— sin Of sincos 77 
apaad - --, (A9) 


sin^ 77 


where the function / = / ( 5 , 6>o) for both types of the 
considered SLs is defined be Eq. (|20l). 
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